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In recent years, much attention has been paid to the study of network structures 
underlying artifactual as well as natural complex systems. There, a vertex and an edge 
represent an element and a relation between elements, respectively. A main interest here 
is in the topological structure of the network, and the physical details of interactions 
between elements play rather marginal roles. This certainly offers an efficient approach 
to understanding collective behaviors of complex systems. The study along this line was 
initiated by the works of Watts and Strogatz [1] on small- world networks and of 
Barabasi and Albert [2] on scale-free networks (see also [3]). These models are 
essentially different from random graphs considered by Erdos and Renyi [4] (see also 
[5]) least in the following two points: i) a small-world network has a value of the 
clustering coefficient [1], which is much larger than that of a random graph, showing 
strong correlation between the edges, and ii) the degree distribution of a scale-free 
network has the form [2] 

with vertex connectivity, k, and constants, A;oe(0, 1) and 7>1, whereas it is the 
Poisson distribution for a random graph. 

The scale-free degree distribution in Eq. (1) implies that there exist a significant 
number of vertices that have high values of connectivity (i.e., hubs). Such a structure is 
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profoundly relevant to the concept of robustness and vulnerability of the network [6]. 

In [2], it has been discussed that it is sufficient for realizing a scale-free network to 
consider its growth with the so-called preferential attachment rule, that is, a vertex 
created anew tends to be linked to an old vertex with probability proportional to 
connectivity of that old vertex. This may be able to explain the structures of the Internet 
[7], the WWW [3], and pattern of citing papers [8]. It, however, does not seem to 
explain the natures of biological networks such as food webs [9], metabolic networks 
[10], and protein architectures in cells [11], for example. This implies that there should 
be a variety of mechanisms, which can lead to scale-free statistics in Eq. (1). 

There are in fact some methods of generating scale-free networks, which do not 
assume growth and preferential attachment. For example, in a static model given in [12], 
a weight factor of the specific form is assigned to each vertex. Then, two vertices are 
randomly selected with the probability proportional to their weights and are connected 
by an edge if they were not already linked (otherwise, they are discarded). Such an 
algorithm leads to Eq. (1) with 7 > 2 for large k. Another static model proposed in [13] 
introduces varying vertex fitness. This idea has further been elaborated in [14] (see also 
[15]). Nowadays, varying fitness is referred to as the hidden variable in the literature. 

In this paper, we construct static complex networks by making use of fluctuating 
random graphs. The discussion may be somewhat similar to the works in [12-15], but it 
is in the spirit of superstatistics [16-20] that establishes a connection between 
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Boltzmann-Gibbs statistical mechanics and its generalizations based on temperature 
fluctuations. An essential point here is the analogy of Erdos-Renyi random graph theory 
to Boltzmann-Gibbs statistical mechanics: in the latter, temperature of a canonical 
ensemble is fixed, while the probability of connecting each two vertices is fixed in the 
former. Therefore, temperature fluctuation corresponds to varying probability of 
connection. Applying the quantum mechanical method, we derive a general formula for 
the hidden variable distribution, which relates a random graph to a general complex 
network. As an example, we analyze a scale-free network, in detail, and present the 
explicit formula for the associated hidden variable distribution. Such a hidden variable 
distribution is shown to decay as a power law. Thus, the numerical result given in [13] 
is confirmed analytically. 

Let us recall the classical discussion about random graphs of Erdos and Renyi. 
Consider A'^ vertices. Each pair of vertices is connected by an edge with probability, q. 
New edges are attached in this way, up to the total number, A'^ - 1 . Then, the probability 
of finding a vertex with degree (i.e., connectivity) k is given by the binomial distribution. 
In the large limit, such a degree distribution becomes Poissonian: 



\^dom(k) = — e-\ (2) 



where A is a fixed constant, q(N-l). This is the analog of the canonical distribution 
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in Boltzmann-Gibbs statistical mechanics, and thus A corresponds to temperature. 

Now, a question is how the degree distribution can change if A is a stochastic 
variable, in analogy with temperature fluctuation in superstatistics [16-20]. In this case, 
the distribution in Eq. (2) is regarded as a conditional probability distribution, 
Prandom I Thus, wc havc 



where p(k) is a certain degree distribution and IT (A) is the associated hidden variable 
distribution [13-15]. In other words, p(k) is the marginal of the joint probability, 
n(A) PrandomC^ I '^)- ^q. (3) Implles that p(k) is the Poisson transformation of the 
hidden variable distribution. In this way, we are able to construct a wide class of 
complex networks from a fluctuating random graph. 

First, let us consider a general form of p(k). This is of obvious importance since 
real-world networks have diverse structures, but they can still be characterized by their 
degree distributions [21]. To find the form of 11 (A), we employ the quantum 
mechanical method. In particular, we use the coherent state of the harmonic oscillator 
with unit mass and frequency. The Planck constant is also set equal to unity. The 
creation and annihilation operators, and d, satisfy the algebra: [d,d^] = l, 
[d,a\ = [d\d^] = 0. The ground state, |0), is annihilated by d: a|0) = 0. The k- 




(3) 
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particle state is given by \k) = (kl) ''^(a' )*|0), which is the eigenstate of the number 
operator, n = a^a, that is, n\k) = k\k). {\k)} ^ ^ forms an orthonormal complete set 

termed the Fock basis. The coherent state, \a), is defined by, 
|a) = e~ "' '^^j^_o(^* / |^)> where a is a complex variable. It also forms the 

(over)complete set, satisfying ;r || | a){a | = 1 , where d^a = d (Re a)d (Im a) and 

the domain of integration is the whole complex a -plane. It should be noticed that the 

I 1 2 /i 1 2 \ ^ 

number distribution in the coherent state is Poissonian: |(/^|«)| = (|a| 1 I k\ 
Now, we consider the P-representation of a density matrix p [22,23]: 
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p = jjj-aP(a, a*)|a){a|. (4) 



The function, P{a, a*), may be singular and negative, in general (however, as we shall 
see, it turns out to behave well as a probability distribution in the present context). 
Assume that p is a certain function of the number operator. 



p = p{n). (5) 



Then, from Eq. (4), it follows that 
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^^d^aP(a, a*) 



nr.- 

A:! 



= p(k). 



(6) 



In this representation, the degree distribution is the analog of the energy (i.e., particle 
number) distribution. It would be of interest to examine possibility of applying the 
second quantization formalism to constructing a network based on such an analogy. 
Since p is the function of the number operator, P(a, a*) is a function only of \cc\^ . 
Therefore, we see the parallelism between Eqs. (3) and (6) with identification, 
A <-> I a I . A crucial point is that it is possible to express P(a, a*) in Eq. (4) 
conversely in terms of p as follows [23]: 



From this equation, we can arrive at the general analytic formula for the hidden variable 
distribution 



P(a, a*) = 



1 



al3*-a*l3 



(7) 



n(X) = 7i;P(a, a*) 




(8) 



with A, = \a 
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Finally, let us calculate as an important example the hidden variable distribution for a 
scale-free network characterized by the degree distribution in Eq. (1). In this case, the 
normalized density matrix is taken to be 



(9) 
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where ^(s, a) is Hurwitz's generalized zeta function [24], which is well defined for 
s>l (but can be analytically continued to an arbitrary complex s except for the 
singularity at 5 = 1). Using the Mellin transformation of p^j, that is, 

Psf ~ [C(7' ^o) r(7)] ^^dtt^'^e with Euler's gamma function T{z), we 

obtain the following analytic formula for the hidden variable distribution: 

The first and second moments of this distribution are analytically calculated to be 



"0' 



(12) 




C(r-2,fco) 



-(2^0 + 1) 



ay-i.fco) 



+ ^o(^o + l)' 



(13) 



respectively. From Eq. (11) as well as these, it is clear that n^f(A) decay as a power 
law 



In Fig. 1, we present the plot of 11^,^(2) with respect to A. There, one recognizes the 
asymptotic power-law behavior of Yi^^{X). This, in turn, explains why the authors of 
[13] could numerically realize a scale-free network by using a power-law hidden 
variable distribution. 

In conclusion, we have constructed in the spirit of superstatistics a generic complex 
network by introducing fluctuation to a random graph. We have derived the exact 
analytic formula for the hidden variable distribution, which describes the fluctuation 
and generates a general form of the degree distribution through the Poisson 



(14) 



with the exponent 



5 = Y. 



(15) 
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transfonnation. As an example, we have considered a static scale-free network. We have 
explicitly calculated the associated hidden variable distribution and have shown that it 
decays as a power law with the same exponent as that of the degree distribution. 

S. A. was supported in part by the Grant-in- Aid for Scientific Research of Japan 
Society for the Promotion of Science. S. T. acknowledges support from the FWF Grant 
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Figure caption 

Fig. 1 The log-log plot of 11^^ (A) with respect to A with 5{=Y) = 2.5 and 
= 0.5. All quantities are dimensionless. 
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